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ABSTRACT  hs 

i 

mi  improved  convection  heat-transfer  model  has 
beer,  developed  for  the  prediction  of  the  transient  ht 

tooth  surface  temperature  of  spur  gears.  The  dissipa¬ 
tive  quality  of  the  lubricating  fluid  is  shown  to  be 
limited  to  the  capacity  extent  of  the  thermal  boundary  J 

layer.  This  phenomenon  can  be  of  significance  in  the 
determination  of  the  thermal  limit  of  gears  accelerat-  k 

ir.g  to  the  point  where  gear  scoring  occurs.  Steady- 
state  temperature  prediction  is  improved  considerably  N 

through  the  use  of  a  variable  integration  time  step 
that  substantially  reduces  computer  time.  Computer-  n 

generated  plots  of  temperature  contours  enable  the 
user  to  animate  the  propagation  of  the  thermal  wave  as  Pd 

the  gears  come  into  and  out  of  contact,  thus  contribut¬ 
ing  to  better  understanding  of  this  complex  problem.  q 

This  model  has  a  much  better  capability  at  predicting 
gear-tooth  temperatures  than  previous  models.  **o 

NOMENCLATURE  *«'  r 


heat-transfer  coefficient  for  gear  sides, 
W/hr  m-  K  (Btu/hr  ft2  °F) 

heat-transfer  coefficient  for  unlubricated 
flank  of  gear*  W/hr  m-  K  (Btu/hr  ft^  °F) 

heat-conversion  factor 

thermal  conduc t ivi ty ,  W/m  K  (2tu/ft  °F) 
number  of  teeth 

normal  coordinate  to  tooth  profile*  m  (in.) 

diameter  to  pitch  ratio,  m~^  (in.“M 

heat  flux,  W/hr  (Btu/hr) 

addendum  radius  of  gear*  m  (in.) 

radius  vector  of  fluid  particle,  m  (in.) 


A  heat  flux  area,  m^  (in.^) 

a  acceleration  of  fluid  particle,  m/sec^ 

(  in . / sec  ^  ) 

b  Herzian  contact  width,  m  (in.) 

C  heat  capacitance  matrix 

Ga  dimensionless  constant  (2R0/F)’,  see  Eq.  (9) 

F  face  width  of  tooth  in  contact,  m  (in.) 

f  coefficient  of  friction 

G  centrifugal  acceleration,  m/sec^  (in./sec^) 

hi  heat-transfer  coefficient  for  lubricated  flank 

of  gear,  W/hr  m^  K  (Btu/hr  ft^  ° F ) 

hea t -transfer  coefficient  for  top  land  of 
gear,  W/hr  m^  K  (Btu/hr  ft^  ®F) 


S  surface  finish 

s  tangential  coordinate  to  the  tooth  profile,  m 

( in .  ) 

T  temperature,  K  (°F) 

t  time*  hr 

Vs  instantaneous  sliding  velocity,  m/s ec  (ft/sec) 

Vt  instantaneous  total  velocity,  m/sec  (ft/sec) 

W  normal  tooth  load  intensity,  N/m  (lb/ft) 

x,y,z  dimensionless  Cartesian  coordinates  of  gear 

tooth 

ct  thermal  diffusivity,  m *Vs  (in.-/s) 

y  finite  difference  weighting  factor  (0  <  y  <  1) 


1 


e 


dimensionless  gear  temperature 
p  dynamic  viscosity 

v  kinematic  viscosity 

co  angular  velocity  vector,  sec’^- 


Considering  a  frame  of  reference,  attached  to  the 
gear  at  its  center  of  rotation  and  thus  rotating  with 
an  angular  velocity  w  (assumed  constant),  the  accel¬ 
eration  of  the  fluid  particle  is  given  by 

a  =  -toxcoxr-2uxV  (1) 


INTRODUCTION 

Recent  trends  in  the  design  of  high-speed  machin¬ 
ery  necessitated  careful  consideration  of  factors 
affecting  gear  durability.  As  a  major  failure  mode  of 
high-speed  gearing,  scoring  is  receiving  significant 
attention.  Scoring  failures  have  been  attributed  to  a 
sudden  surge  of  the  contact  temperature  of  the  gear. 
This  contact  temperature  is  directly  related  to  the 
bulk  temperature  of  the  gear  blank.  Thus,  an  accurate 
prediction  of  scoring  thresholds  hinges  on  a  correct 
evaluation  of  the  bulk  temperature.  To  evaluate  the 
contact  temperature,  finite-element  programs  were 
developed  by  Patir  (1977),  Townsend  (19SI),  and 
El-Bayoumy  (1984)  with  the  purpose  of  analyzing  the 
steady-state  thermal  behavior  of  spur  gears.  The 
steady-state  approximation  required  the  assumption  of 
constant  heat-transfer  coefficients.  The  numerical 
values  of  the  coefficients  were  based  on  correlation 
with  experimental  measurements  made  by  Townsend  (1981). 
When  operating  conditions  sue  i  as  speed,  applied  load, 
initial  temperature,  and  oil  _et  velocity,  orossure, 
and  viscosity  were  varied,  correlation  between  theory 
and  experiment  was  not  always  *, at  is  factory . 

The  difficulties  arising  from  havin?  to  assume  the 
proper  values  for  these  heat-transfer  re*  .f  i:  ic  ients  were 
overcome  in  El-Bayoumy  (1985),  where  a  transient  ther¬ 
mal  analysis  approach  was  developed.  The  convection 
process  in  El-Bayoumy  (1985)  was  based  on  an  extended 
model  from  Van  Heijningen  (1974)  that  yielded  a  set  of 
time-dependent  coefficients. 

Closer  examination  of  the  hydrodynamic  and  thermal 
phenomena  within  the  adjacent  fluid  film  shows  the  film 
thickness  remaining  above  the  impingement  point  to 
decrease  with  time.  The  thermal  boundary-layer  thick¬ 
ness,  on  the  other  hand,  grows  steadily  with  time.  An 
instant  is  therefore  reached  when  the  two  thicknesses 
are  equal.  At  that  point  the  fluid  film  is  fully  satu¬ 
rated  with  heat  and  can  no  longer  participate  in  the 
convection  process.  The  wetted  side  of  the  gear  tooth 
may  therefore  be  divided  into  two  regions:  a  convec¬ 
tive  region,  where  the  heat-transfer  coefficients, 
derived  from  the  extended  Blok  model,  can  be  used  and 
a  nonconvec t ive  region  with  vanishing  heat-transfer 
coef  f icients . 

The  objective  of  this  paper  is  to  quantify  the 
influence  of  the  interaction  between  the  thermal  bound¬ 
ary  and  the  oil  film  thickness  on  the  transient  thermaL 
behavior  of  spur  gears.  In  addition,  the  thermal  wave 
propagation  process  of  the  hot  spots  are  portrayed  by 
instantaneous  isotherms  as  they  move  through  the  gear 
tooth.  The  asymptotic  approach  towards  steady  state  in 
the  absence  and  presence  of  the  thermal  boundary  layer 
is  evaluated. 

Interaction  Between  Thermal  and  Hydrodynamic  Effects 
Within  the  Fluid  Film 

To  evaluate  how  much  of  the  adjacent  lubricant 
actually  contributes  to  the  heat  dissipation  of  the 
gear  tooth,  :r.e  need.*  l o  consider  the  hydro**/ namic 
forces  acting  on  a  fluid  particle  near  the  tooth  sur¬ 
face  above  and  below  the  point  of  impingement.^- 


^•The  derivations  given  in  this  section  are  based  on 
private  communications  with  Dr.  John  Murdock  of  the 
Aerospace  Corporation. 


where  r  is  the  radius  vector  to  the  particle,  and 
V  is  the  relative  velocity  of  jet  impingement.  The 
two  terms  appearing  on  the  right  h^nd  side  of  Eq.  (1) 
are  the  centrifugal  and  coriolis  acceleration  compo¬ 
nents.  Since  r  varies  very  little  across  the  tooth, 
the  centrifugal  acceleration,  for  all  practical  pur¬ 
poses  can  be  considered  constant  for  points  along  the 

9 

tooth  profile.  Its  magnitude  is  V“/r  where  V  is 

8  8 

the  velocity  of  the  gear,  at  the  pitch  point  the  cen¬ 
trifugal  acceleration  as  seen  from  the  sign  in  Eq.  (1) 
tends  to  strip  off  the  oil  layer  from  the  surface  of 
the  tooth.  The  magnitude  of  the  second  term,  the 
coriolis  acceleration,  is  2VgVj/r  where  Vj  is  the 
magnitude  of  the  jet  velocity  relative  to  the  gear. 
Since  Vj  is  of  the  same  order  of  magnitude  as  V, , 
the  coriolis  component  cannot  be  neglected.  For 
points  on  the  tooth  profile  below  impingement,  the 
direction  of  tne  coriolis  acceleration  is  normal  to 
the  tooth  surface  pointing  outward.  Conversely,  for 
points  above  impingement  the  coriolis  acceleration 
points  inward,  as  illustrated  in  Fig.  1,  establishing 
a  "pseudo-gravitational"  field.  Thus,  the  initial 
area  A  of  convection  heat  flux  q,  characterized  by 
the  equation 


q  =  hA  AT  (2) 

covers  the  length  of  the  tooth  profile  above  the 
impingement  point.  How  much  of  this  area  will  partici¬ 
pate  in  heat  dissipation  can  be  determined  when  the 
thermal  and  fluid  boundary  layers  are  examined  versus 
time.  The  thicknesses  of  these  layers  as  a  function 
of  time  are  schematically  illustrated  in  Fig.  2.  The 
viscous  layer  develops  in  a  time  of  the  same  order  of 
magnitude  as  the  oil  spread  time  which  is  negligible 
compared  with  the  heat-transfer  time.  While  the  vis¬ 
cous  forces  tend  to  attach  the  oil  film  to  the  tooth 
surface,  centrifugal  acceleration  tries  to  strip  it 
off.  The  amount  of  film  remaining  on  the  surface  has 
been  determined  by  Jefferis  (1967),  who  showed  that  the 
thickness  of  this  oil  film  is  given  by 


H 


f 


(3) 


where  v  is  the  kinematic  viscosity  of  the  oil,  x  is 
the  coordinate  along  the  tooth  profile  measured  from 
the  impingement  point  (see  Fig.  1),  G  is  the  centrifu¬ 
gal  acceleration,  and  t  is  time.  This  relationship 
is  depicted  in  Fig.  2  as  a  straight  line  with  negative 
slope . 

The  thermal  boundary  layer  thickness  is  growing  as 
a  function  of  time  according  to  the  deWinter  and  Bloc 
(1974)  formula 


Ht  =  1.362^/Jt  (4) 

where  a  is  the  thermal  diffusivity  of  the  oil.  The 
distance  from  the  point  of  impingement  of  the 

point  at  which  the  thermal  boundary  layer  thicknesses 
of  Eq.  (3)  is  equal  to  the  film  thickness  of  Eq.  (4) 
and  is  therefore  given  by 


x 


1 


1.85 


Ga  2 


v 


(5) 
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This  is  Che  point  cf  heat  saturation  of  the  fluid  film. 
For  points  on  the  tooth  profile  below  this  point  (i.e., 
x  ^  xi),  the  fluid  no  longer  contributes  to  the  convec¬ 
tive  dissipation  of  heat.  In  other  words,  point 
becomes  the  effective  impingement  point.  For  points 
above  P i. ,  a  heat-transfer  coefficient  based  on  the 
solution  of  the  one-dimensional  heat  equation 


o 


(6) 


is  used.  This  coefficient  may  be  written  as 

hj  =  K/(irat)1/2  (7) 

Mathematical  Formulation  of  the  Transient  Thermal 
Problems  in  Gearing 

As  noted  in  El-Eayoumy  (1985),  the  partial  differ¬ 
ential  equation  governing  heat  transfer  within  the  gear 
sector  shown  in  Figs.  3  and  i  is  given  by 


where 


n  n/R0 

n  normal  coordinate  to  the  tooth  profile 

hj  convective  heat  transfer  coefficient  in  tooth 

profile  above  the  effective  point  of  impinge¬ 
ment,  W/hr-m2  K  (Btu/hr  ft2  *F) 

hs  convective  heat  transfer  coefficient  on  front 
and  back  of  gear 


hjt  convective  heat  transfer  coefficient  on  top  land 

ht  convective  heat  transfer  coefficient  on  remain¬ 

der  of  tooth  including  bottom  land 

Pj  dimentral  pitch  of  gear,  (in.-'-) 


q*  q/qav 


,2  36  n 
aR  r —  =  u 
o  dt 


Tj  tooth  involute  boundary  above  saturation  point 
rs  front  and  back  end  faces  of  tooth 

rc  remainder  of  tooth  profile 


whe  re 


x 


and  where 


>  (9) 


e  = 


k?  (T  -  T  ) 
d _  o 


It  was  possible  in  El-Bayoumy  (1985)  to  reduce  the 
above  boundary  value  problem  through  discretization  to 
the  recurrence  relation 

(C/At  +  Hy)9n+1  =  [C/At  -  H(1  -  2 )  ]©n  +  P  (11) 

where  C  and  H  are  the  capacitance  and  heat  conduc¬ 
tivity  matrices,  6n  is  the  dimensionless  temperature 
vector  at  time  tn,  and  0  <  y  <  1.  As  noted  in 
Zienkiewicz  (1977),  for  numerical  stability  y  >  1/2. 

Frictional  Heat  Vector  Options 

The  instantaneous  heat  generated  during  contact 
between  the  two  gears  may  be  written  as 


T  current  gear  temperature,  ”C  (”F) 

T0  initial  gear  temperature,  °C  (°F) 


V 


s 


(12) 


a  gear  diffusivity,  m2/s  (in.2/s) 

k  thermal  conductivity  of  gear  material,  W/m  K 

(Btu/ft  °F) 

R3  outside  radius  of  gear,  m  (ft) 

F  face  width  of  gear,  m  (ft) 

qav  average  generated  heat  flux  through  the  gear  pro¬ 
file,  W/hr  (Btu/hr) 


where 

J  mechanical  equivalent  of  heat,  W/Nm  (Btu/ftlb) 
W  normal  load  intensity,  N/m  (lb/ft) 
b  width  of  contact  region,  m  (ft) 

f  coefficient  of  friction 

Vs  instantaneous  sliding  velocity,  m/s  (ft/s) 


The  associated  boundary  conditions  are 


h . (s , t )R 


ae  “f3' 

— —  *  — 1 — - 


9  -  P.R  q*  =  0  on  T. 
d  o 


-a  h  (s  ,  t  )R 
o9  t _  o 


3n 


9=0  on  T 


>  (10) 


FQ  hs(t)F 

~  <■  -  9  =  0  on  T 

<jz  If.  s 


Assuming  an  even  split  between  the  pinion  and  the  gear, 
the  frictional  heat  input  to  either  gear  is 


U 


TTZ  fWV 
2Jb  s 


(13) 


Various  options  from  the  literature  (Anderson,  1981) 
are  available  for  the  calculation  of  the  friction  ccef 
ficient.  The  following  choices  are  incorporated  into 
the  program  (according  to  Kelly-Benedict,  1961): 


fL  =  O.Clz/ 


log  / 3.17x10  W 


f^svt 


(li) 


f 


""N 


ith _ L 


cn 


If  the  effect  of  surface  finish  is  to  be  examined,  a 
modified  Kelly-Benedict  formula  is 


The  results  presented  in  this  paper  are  based  on  the 
Keiley-Benedict  formula  (Eq.  (14))  as  it  correlates 
well  with  efficiency  measurements  given  in  Anderson 
(1931). 

RESULTS  AND  DISCUSSIONS 

The  heat-transfer  coefficients  derived  in  Eq.  (7) 
were  used  in  El-3ayoumy  (1935)  for  all  points  above  the 
true  impingement  point.  The  computer  program  developed 
for  the  purpose  of  evaluating  transient  thermal  behav¬ 
ior  was  mainly  concerned  with  numerical  stability  of 
two-  and  three-point  finite-difference  schemes.  The 
computer  program  was  also  developed  for  a  comparative 
study  of  gear  temperatures  basea  on  heat-transfer  coef¬ 
ficients  used  in  the  steady-state  program  with  those 
based  on  Eq.  (7).  Preliminary  results  were  given  in 
the  form  cf  snap  shots  of  gear  temperature  contours  at 
the  initial  stages  of  tooth  engagement.  Eecause  the 
program  does  many  iterations,  it  requires  considerable 
time  to  reach  the  steady-state  equilibrium.  Because 
all  the  experimental  data  are  at  steady-state  condi¬ 
tions,  ir  was  necessary  to  run  the  program  for  some 
time  to  match  the  experimental  data. 

The  case  investigated  here  is  that  of  a  pair  of 
eight-pitch  standard  spur  gears  with  28  teeth  each  run¬ 
ning  at  10  0C0  rpm  with  a  normal  load  per  unit  face 
width  of  5912  N/cm  (3378  lb/in.).  The  actual  jet 
impingement  depth  is  assumed  at  87.5  percent  of  the 
tooth  height.  The  resulting  temperature  history  is 
displayed  in  three  forms.  First,  at  predefined  time 
instants,  a  scan  is  made  over  all  elements  for  the  max¬ 
imum  and  minimum  temperature  rise.  Then  the  range  is 
divided  into  a  user-defined  number  of  equal-temperature 
increments.  Each  division  is  thus  associated  with  a 
temperature  level.  These  levels  can  be  presented  in 
the  form  of  contour  plots  such  as  those  shown  in 
rig.  5.  The  analysis  is  first  carried  out  neglecting 
boundary-layer  interaction  (assuming  that  the  fluid 
film  continues  to  absorb  heat.)  The  jet  is  assumed  to 
impinge  5’  past  the  last  point  of  tooth  engagement. 

The  temperature  range  at  each  instant  is  divided  into 
10  equal  divisions,  or  contour  levels.  Because  the 
gear  temperature  will  continuously  vary  with  time,  each 
contour  level  will  define  a  different  temperature  level 
each  time.  For  example,  the  first  contour  level  defines 
a  temperature  rise  of  26.4  °C  (47.5  °E)  at  0.0167  msec 
Cl’),  35.1  aC  (63.1  *F)  at  0.05  msec  (3°),  and  8.3  *C 
(15.2  *F)  at  0.167  msec  (10“)  after  the  beginning  of 
tooth  engagement. 

In  order  to  better  visualize  the  chronological 
development  of  the  heat  wave,  contour  levels  are 
initially  preassigned,  and  the  progress  of  each  contour 
is  monitored  step  by  step.  Thus  the  same  temperature 
history  may  be  displayed  in  form  of  the  fixed  contour 
levels  shown  in  Fig.  6.  The  range  of  the  temperature 
r;se  during  the  first  cycle  of  tooth  engagement,  having 
been  determined  as  approximately  75  *C  (135  ”F),  is 
divided  into  10  equal  divisions  starting  with  the  first 
contour  level  at  75  *C  (135  *F)  and  ending  with  0  "C 
(0  *F)  for  the  last  contour  level.  The  progress  of 
heat-wave  propagation  can  easily  be  traced  through  the 
sequence  of  the  snap  shots  as  the  gear  tooth  moves 
through  engagement  given  in  Fig.  6.  First,  lower  level 
contours  emerge  then  proceed  through  the  gear  tooth, 
followed  by  the  appearance  of  higher  level  contours 
then  a  gradual  withdrawal  of  these  contours  occurs  as 
the  contact  nears  the  pitch  point.  The  contours  then 


reemerge  as  the  tooth  mesh  proceeds  through  the  angle 
of  recess  until  all  the  10  contours  are  realized  at  the 
end  of  contact.  A  significant  drop  in  temperature  then 
follows  due  to  the  absence  of  frictional  heat  and  the 
continued  dissipation  of  heat  through  the  fluid  film. 

The  results  may  be  summarized  in  a  single  tempera¬ 
ture  rise  versus  time  graph  (shown  in  Fig.  7)  illus¬ 
trating  how  the  maximum,  average  midplane,  and  average 
face  temperatures  on  the  loaded  and  coast  sides  of  the 
tooth  vary  with  time.  As  indicated  in  the  figure,  most 
of  the  thermal  activity  is  limited  to  the  duration  of 
tooth  engagement.  Furthermore  the  highest  temperature 
is  highly  localized  at  the  tip  of  the  driving  gear  and 
at  a  point  between  the  lowest  point  of  contact  and  the 
pitch  point.  In  the  beginning  of  tooth  engagement  the 
tip  is  significantly  hotter  than  the  other  localized 
hot  spot. 

The  influence  of  boundary-layer  interaction,  as 
illustrated  in  Figs.  8  to  10,  accentuates  that  differ¬ 
ence.  This  can  be  expected,  because  at  the  lower  hot 
spot  the  fluid  film  continues  to  convect  heat  as  the 
thermal  boundary  layer  has  not  fully  developed  as  it 
will  when  engagement  approaches  the  tip.  Thus,  while 
the  lower  hot  spot  goes  through  a  temperature  rise  of 
40.75  °C  (73.35  *F)  during  initial  engagement,  the  tip 
can  rise  by  as  much  as  96.02  °C  (172.84  °F)  (Figs.  8 

and  9)  versus  77.24  °C  (139.02  *F)  (Figs.  5  and  6)  if 
boundary  layer  interaction  is  neglected.  Note  that  the 
lower  hot  spot  is  unaffected  by  boundary-layer  interac¬ 
tion  (compare  Figs.  5(a)  and  8(a)). 

When  the  simulation  process  is  continued  to  steady 
state,  significant  differences  in  predicted  tooth  tem¬ 
peratures  can  arise  when  boundary-layer  interaction  is 
not  neglected.  This  is  clearly  illustrated  in  Figs.  11 
to  16.  The  maximum  temperature  rise  predicted  with 
boundary-layer  interaction  is  119  °C  (215  *F)  (Figs.  14 
and  15)  in  close  agreement  with  experimental  findings 
shown  in  Fig.  12  from  Townsend  (1981),  compared  with 
78  °C  (141  °F)  (Figs.  11  and  13)  without  boundary-layer 
interaction.  The  errors  are  more  significant  when 
average  tooth  temperatures  are  compared.  With  boundary- 
layer  interaction  the  average  steady-state  temperature 
rises  49  "C  (88  *F)  (Fig.  15).  Without  boundary-layer 
interaction  a  predicted  average  steady-state  tempera¬ 
ture  rise  of  only  14  »c  (25  *F)  is  achieved  (Fig.  13). 

In  Fig.  16  the  maximum  temperature  rise  is  plotted 
versus  time  over  the  entire  simulation  period.  It  is 
interesting  to  note  that  with  boundary— layer  effects 
the  temperature  attains  a  maximum  before  steady  state 
is  reached  and  scoring  could  take  place  even  though  the 
indicated  steady-state  temperature  is  below  the  scoring 
threshold  of  the  tooth  pair. 


SUMMARY  OF  RESULTS 

The  interaction  between  the  thermal  and  lubricant 
boundary  layers  has  been  incorporated  into  a  transient 
thermal  analysis  program  for  predicting  surface  temper¬ 
atures  in  high-performance  spur  gears.  As  a  result  of 
this  interaction  a  portion  of  the  fluid  film  adjacent 
to  the  tooth  becomes  saturated  with  heat.  This  portion 
continues  to  spread  with  time  until  it  covers  the  entire 
wetted  area.  Predicted  temperature  rise  for  a  given  set 
of  gears  was  found  to  correlate  well  with  previously 
measured  data.  If  boundary  layer  interaction  is 
ignored,  the  predicted  temperature  rise  at  the  hottest 
spot  could  be  off  by  as  much  as  42.5  percent.  Two  hot 
spots  exist  on  the  tooth  surface  of  the  driving  gear: 
at  the  tip  of  the  loaded  side  and  below  the  pitch  point. 
The  gear  reaches  its  hottest  temperature  before  steady- 
state  conditions  are  achieved.  Through  the  use  of  pre¬ 
scribed  temperature  contours,  the  process  of  heat-wave 
propagation  through  the  gear  tooth  can  be  animated. 
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FIGURE  1.  -  CORIOLIS  ACCELERATION  ABOVE  AND  BELOW  IMPINGEMENT 
ON  GEAR  TOOTH. 
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FIGURE  2.  -  CHARACTERISTIC  LENGTH  AND  TIME  SCALES  FOR  OIL 
FILM  ON  ROTATING  GEAR  TOOTH. 
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FIGURE  3.  -  CONVECTIVE  HEAT 
TRANSFER  COEFFICIENT  ZONES. 


FIGURE  A.  -  FINITE  ELEMENT 
GRIO. 
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BTTWfIN  MINIMUM  ANII  MAXIMUM  UMPERATURE  .  SPITI)  10  000  rpm; 
l  OAI)  Oil?  N/cm  (33/8  Iii/IN.T;  8/. 6  PFRCFN1  OIL  JET  IMPINGEMENT 
DEPTH;  NO  BOUNDARY  LAYER  EFFECT . 
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<e)  in0-  (O  18°. 


(g)  M6.50* 

FIGURE  6.  -  TEMPERATURE  RISE  CONTOURS  (IN  DEGREES  CELSIUS)  AT  1° 
TO  96.5°  OF  GEAR  ROTATION  AFTER  INITIAL  TOOTH  CONTACT.  INITIAL 
NONEQUILIBRIUM  CONDITIONS  WITH  10  EQUAL  PREASSIGNED  FIXED  CON¬ 
TOUR  LEVELS  BETWEEN  MINIMUM  AND  MAXIMUM  TEMPERATURE.  SPEED 
10  000  RP*;  LOAD  5912  N/cm  (5578  LB/IN.):  87.5  PERCENT  OIL  JET 
IMPINGEMENT  DEPTH;  NO  BOUNDARY  LAYER  EFFECT. 
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FIGURE  15.  -  TEMPERATURE  RISE  VERSUS  TIME  OR  GEAR  ROTATION  AFTER  65  REVOLUTIONS. 

SPEED.  10  000  RPM;  LOAD.  5912  N/cm  <5578  LB/IN.);  87.5  PERCENT  OIL  JET  IMPINGE¬ 
MENT  DEPTH;  WITHOUT  BOUNDARY  LAYER  EFFECTS  AFTER  STEADY-STATE  EQUILIBRIUM  CONDITIONS. 
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FIGURE  15.  -  TEMPERATURE  RISE  VERSUS  TIME  OR  DEGREES  ROTATION  AFTER  %  REVOLUTIONS.  SPEED. 
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FIGURE  16.  -  PROFILE  MAXIMUM  TEMPERATURE  VERSUS  GEAR  REVOLUTIONS  WITH  AND  WITHOUT  BOUNDARY 
LAYER  EFFECTS. 
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